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A critical discussion of the present signals for the phase transition to quark-gluon plasma 
(QGP) is given. Since hadronic rescattering models predict much larger flow than ob- 
served from 1 to 50 A GeV laboratory bombarding energies, this observation is inter- 
preted as potential evidence for a first-order phase transition at high baryon density. A 
detailed discussion of the collective flow as a barometer for the equation of state (EoS) of 
hot dense matter at RHIC follows. Here, hadronic rescattering models can explain < 30% 
of the observed elliptic flow V2 for pT > 2 GeV/c. This is interpreted as an evidence for 
the production of superdense matter at RHIC. The connection of V2 to jet suppression is 
examined. A study of Mach shocks generated by fast partonic jets propagating through 
the QGP is given. The main goal is to take into account different types of collective 
motion during the formation and evolution of this matter. A significant deformation of 
Mach shocks in central Au+Au collisions at RHIC and LHC energies as compared to the 
case of jet propagation in a static medium is predicted. A new hydrodynamical study of 
jet energy loss is presented. 

1. Observables for the QGP phase transition 

Lattice QCD calculations yield a phase diagra (see Fig. [Ij which shows a 
crossing, but no first-order phase transition to the quark-gluon plasma (QGP) for 
vanishing or small chemical potentials /is, i.e. for conditions accessible at central 
rapidities at full RHIC energy. A first-order phase transition is expected to occur 
only at high baryochemical potentials or densities, i.e. at the lower SPS and RHIC 
energies and in the fragmentation regi on o f RHIC, y w 3 — 5 . Here, the criti- 
cal baryochemical potential is predictecfDl to be large , fi% « 400 ± 50 MeV, and 
the critical temperature to be Tc ~ 150 — 160 MeV. We expect a first-oder phase 
transition also at finite strangenesJ^. Predictions for the phase diagram of strongly 
interacting matter for realistic non- vanishing net strangeness are urgently needed to 
obtain a comprehensive picture of the QCD phase structure in all relevant dimen- 
sions (isospin, strangeness, non-equilibrium) of the EoS. Multi-strange degrees of 
freedom are very promising probes for the properties of the dense and hot matteJ^. 
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Fig. 1. The phase diagram with the critical end point at ~ 400 MeV, T 160 MeV as predicted 
by Lattice QCD. In addition, the time evolution in the T — ^g— plane of a central cell in UrQMD 
calculations [from Bravina et a/.jSI is depicted for different bombarding energies. Note that the 
calculations indicate that bombarding energies Eji^^b ^ 40 A GeV are needed to probe a first- 
order phase transition. At RHIC this point is accessible in the fragmentation region only [from 
Bratkovskaya et ai.f3. 



1.1. Thermodynamics in the T - plane 

Figure [T] shows a comparison of the QCD predictions with the thermodynamic pa- 
rameters T and hb extracted from the UrQMD transport model in the central 
overlap regime of Au-|-Au collision!^. Full dots with errorbars denote the 'experi- 
mental' chemical freeze-out parameters - determined from fits to the experimental 
yields - taken from Ref.^. Triangular and quadratic symbols (time-ordered in verti- 
cal sequence) stand for temperatures T and chemical potentials /is extracted from 
UrQMD transport calculations in central Au-|-Au (Pb+Pb) collisions at RHIC'^ as 
a function of the reaction time (separated by 1 fm/c steps from top to bottom). 
Open symbols denote nonequilibrium configurations and correspond to T parame- 
ters extracted from the transverse momentum distributions, whereas the full sym- 
bols denote configurations in approximate pressure equilibrium in longitudinal and 
transverse direction. 

During the nonequilibrium phase (open symbols) the transport calculations 
show much higher temperatures (or energy densities) than the 'experimental' chem- 
ical freeze-out configurations at all bombarding energies (> 11 A GeV). These 
numbers are also higher than the critical point (circle) of (2+1) flavor lattice QCD 
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calculations by the Bielefeld-Swansea-collaboratiorP (large open circle) and by the 
Wuppertal-Budapest-coUaboratiorPl (the star shows earlier results frorrP^ . The en- 
ergy density at /Xc, Tc is of the order of w 1 GeV/fm^. At RHIC energies a cross-over 
is expected at midrapidity, when the temperature drops during the expansion phase 
of the 'hot fireball'. The baryon chemical potential /i^ has been obtained from a 
statistical model analysis by the BRAHMS collaboration based on measured an- 
tihadron to hadron ratios^ for different rapidity intervals at RHIC energies. At 
midrapidity one finds /x_b — 0, whereas at forward rapidities /is increases up to 
fiB — 130 MeV at y = 3. Thus only a forward rapidity measurement (y « 4 — 5) 
will allow to probe large at RHIC. The STAR and PHENIX detectors at RHIC 
offer a unique opportunity to reach higher chemical potentials and the first-order 
phase transition region at midrapidity in the HiMu-RHIC-running at = 4—12 
GeV in the coming year. The International FAIR Facility at GSI will be offering a 
fully devoted research program in the next decade. 



1.2. Hydrodynamic flow 

Hydrodynamic flow and shock formation has been proposed 

earljUHni 

as the 

key mec han ism for the creation of hot and dense matter in relativistic heavy-ion 
collisionJISl xhe full three-dimensional hydrodynamical flow problem is much more 
complicated than the one-dimensional Landau model "'^^r the 3-dimensional com- 
pression and expansion dynamics yields complex triple differential cross-sections 
which provide quite accurate spectroscopic handles on the EoS. The bounce-off, 
vi{pt) (i.e., the strength of the directed flow in the reaction plane), the squeeze- 
out, V2{pt) (the strength of the seco nd moment of th e azirnuthal pa rticle emis- 
sion distribution:>1111^ 15 | l« | 17 | 18 | 19| 

and the 

antifloJ15|16|17|18119] 

(third flow 

componenlPQEIl) serve as differential barometers for the properti es of cornpressed, 
dense matter from SIS to RHIC. In particular, it has been 

3howill^|15|16|17|18|19| 

that the disappearance or " collapse" of flow is a direct result of a first-order phase 
transition. 

|00| 

Several hydrodynamic models^ have been used in the past, starting with the 
one-fluid ideal hydrodynamic approach. It is well known that the latter model pre- 
dicts far too large flow effects. To obtain a better description of the dynamics, 
viscous fluid models have been developecf^ ■^^ '^'^ In parallel, so-called 

three-fluid models, which distinguish between projectile, target and the fireball 
fluid, have been considered^Ql. Here viscosity effects appear only between the differ- 
ent fluids, but not inside the individual fluids. The aim is to have at our di sposal a 
reliable, three-dimensional, relativistic three-fluid model including viscositjE^ESl 

Flow can be described very elegantly in hydrodynamics. However, also con- 
sider microscopic multicomponent (pre-)hadron transport theory, e.g. models like 
qMEv31, iQMD^^^ UrQMD or HSE^^^ control models for viscous hydro- 
dynamics and as background models to subtract interesting non-hadronic effects 
from data. If hydrodynamics with and without quark matter EoS, hadronic trans- 
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Fig. 2. Time evolution of directed flow Px/N as a function of rapidity for Au+Au collisions at 
8 A GeV in the one-fluid model. Left: Hadronic EoS without phase tra nsition. Right; An EoS 
including a flrst-order phase transition to the QGP [from BrachmannjSSl 



port models without quark matter - but with strings - are compared to data, can 
we learn whether quark matter has been formed? What degree of equilibration 
has been reached? What does the EoS look like? How are the particle properties, 
self-energies, cross sections changed? 



1.3. Review of AGS and SPS results 

Microscopic (pre-)hadronic transport models describe the fo rma tion and distribu- 
tions of many hadronic particles at AGS and SPS rather welPS. Furthermore, the 
nuclear EoS has been extracted by comparing to flow data which are described 
reasonably well up to AGS energies ■J" ■J' "J" '^-J Ideal hydrodynamical calcu- 

lOOl 

lations, on the other hand, predict far too much flow at these energies^. Thus, 
viscosity effects have to be taken into account. 

In particular, ideal hydrodynamical calculations yield factors of two higher for the 
sideward flow at SIS^"^ and AGS, while the directed flow Px/m measurement of the 
E895 collaboration shows that the p and A data are reproduced reasonably welPSEH 
in UrQMD, i.e., in a hadronic transport theory with reasonable cross-sections, i.e. 



1, 2008 16:46 WSPC/INSTRUCTION FILE paperl2 



Hydrodynamic Flow and Mach Cones 5 




Fig. 3. Measured SIS and AGS proton dpx/dy-slope data compared to a one-fluid hydrodynamical 
calculation. A lin ear extrapolation of the AGS data indicates a collapse of flow at i?£a6 ~ 30 A 
GeV [see also Ref.H^j^ i.e. for the lowest SPS- and the upper FAIR- energies at GSI [from Paech 
et al.pM The point at 40 A GeV is calculated using the NA49 central data, cf. Fig. |4] 



realistic mean-free-path of the constituents. 

Only ideal hydrodynamical calculations predic t, how ever, the appearance of 
a so-called "third flow component'^l^ or " antiflow'SIISll in central collisions. We 
stress that this only holds if the matter undergoes a first-order phase transition 
to the QGP. The signal is that around midrapidity the directed flow, Px{y), of 
protons develops a negative slope! In contrast, a hadronic EoS without QGP phase 
transition does not yield such an exotic "antiflow" (negative slope) wiggle in the 
proton flow vi{y). The ideal hydrodynamic time evolution of the directed flow, 
Px/N, for the purely hadronic EoS (Fig. [2]l.h.s.) does show a clean linear increase 
of Px{y), just as the microscopic transport theory and as the datePSl. For an EoS 
including a first-order phase transition to the QGP (Fig. [2]r.h.s.) it can be seen, 
however, that the proton flow vi ~ Px/pt collapses; the collapse occurs around 
midrapidity. This observation is explained by an antiflow com ponent of protons, 
which develops when the expansion from the plasma sets irl^. 

The ideal hydrodynamic directed proton flow px (Fig. [3|) shows even negative 
values between 8 and 20 A GeV. An increase back to positive flow is predicted 
with increasing energy, when the compressed QGP phase is probed. But, where 
is the predicted minimum of the proton flow in the data? The hydrodynamical 
calculations suggest this "softest-point collapse" is at E^ab ~ 8 A GeV. This has 
not been verifled by the AGS data! However, a linear extrapolation of the AGS 
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Fig. 4. VI at SPS, 40 A GeV and 158 A GeV [from Alt et ai.pSl. The proton antiflow is observed 
in the NA49 experiment even at near central collisions, which is in contrast to the UrQMD model 
involving no phase transition (Fig[5]l- 
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Fig. 5. Proton and pion flow vi = px/pT at 40 A GeV as obtained within the UrQMD model. No 
pro ton antiflow is generated in this hadronic transport theory without phase transition (c.f. Ref. 

m 



data indicates a collapse of the directed proton flow at E^ab ~ 30 A GeV (Fig. [3]). 

Recently, substantial support for this prediction has be en obtained by the low 
energy 40 A GeV SPS data of the NA49 collaboratior^SHSl (pf Fig.|l|). These data 
clearly show the first proton " antiflow" around mid-rapidity, in contrast to the AGS 
data as well as to the UrQMD calculations involving no phase transition (Fig. \^ . 



1, 2008 16:46 WSPC/INSTRUCTION FILE paperl2 



Hydrodynamic Flow and Mach Cones 7 



0.06 
, 0.05 
0.04 
0.03 
0.02 
0.01 





;l , 40 GeV/A 



p, < 2 GeV/c 



158 GeV/A 



50. 



ooo 



I 

' p, < 2 GeV/c 




Fi g. 6 . Elliptic flow V2 of protons versus rapidity at 40 A GeV Pb+Pb collisions [from Alt et 
aZ.p^ as measured by NA49 for three centrality bins: central (dots), mid-central (squares) and 
peripheral (triangles). 



Thus, at bombarding energies of 30-40 A GeV, a first-order phase transition to 
the baryon-rich QGP is most Ukely observed; hence the first order phase transition 
line is crossed (cf. Fig. [1]). This is the energy region where the new FAIR facihty 
at GSI will operate. There are good prospects that the baryon flow collapses and 
other first-order QGP phase transition signals can be studied soon at the lowest 
SPS energies as well as at the RHIC planned HiMu-run at midrapidity as well as 
the fragmentation region y > 4 — 5 for the highest RHIC and LHC-coUider energies. 
These experiments will enable a detailed study of the first-order phase transition 
at high fiB and of the properties of the baryon-rich QGP in the near future. 

2. Proton elliptic flow collapse at 40 A GeV - evidence for a 
first-order phase transition at highest net baryon densities 

At SIS energies, microscopic transport models reproduce the data on the excita- 
tioi i funct ion of the proton elliptic flow V2 quite well: A soft, momentum-dependent 
seems to account for the data. The observed proton flow V2 below ~ 5 
A GeV is smaller than zero, which correspo nds to the squeeze-out predicted by 
hydrodynamics long agcl ^^ l ^^ l l^ l ^6 | 17 | 18 119j 'pj^g ^^^^ exhibit a transition 

from squeeze-out to in-plane flow in the midrapidity region. The change in sign 
of the proton V2 at 4-5 A GeV is in accord with transport calculations (UrQMD 
calculations^^ for HSD results see Ref.'^^ ^^). At higher energies, 10-160 A GeV, a 
smooth increase of the flow V2 is predicted from the hadronic transport simulations. 
In fact, the 158 A GeV data of the NA49 collaboration suggest that this smooth 
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increase proceeds between AGS and SPS as predicted. Accordingly, UrQMD cal- 
culations without phase transition gi ve cons iderable 3% f 2 flow for midcentral and 
peripheral protons at 40 A GeV (RefSIEHl). 

This is in strong contrast to recent NA49 data at 40 A 

Ge\M5J46| (cf. Fig. [6]): A 

sudden collapse of the proton flow is observed for midcentral as well as for peripheral 
protons. This collapse of V2 for protons around midrapidity at 40 A GeV is very 
pronounced while it is not observed at 158 A GeV. 

A dramatic collapse of the flow vi is also observed by again around 

40 A GeV, where the collapse of V2 has been observed. This is the highest energy 
- according to Ref.-'^ and Fig. [1]- at which a flrst-order phase transition can be 
reached at the central rapidities of relativistic heavy-ion collisions. We therefore 
conclude that a first-order phase transition at the highest baryon densities accessi- 
ble in nature has been seen at these energies in Pb-|-Pb collisions. Moreover, Ref.S9l 
shows that the elliptic flow clearly distinguishes between a first-order phase transi- 
tion and a crossover. 



3. Mach shocks induced by partonic jets in expanding QGP 

Sideward peaks have been recently 

observecPEHaSS] 

in azimuthal distributions 

of secondaries associated with the high-px hadrons in central Au+Au collisions at 
^/s = 200 GeV. In Ref.^^ such peaks had been predicted as a signature of Mach 
shocks created by partonic jets propagating through a QGP formed in heavy-ion 
collisi ons. Analogous Mach shock waves were studie d prev iously in cold hadronic 

niatteiiniiiiiisassi 

as well as in n uclear Fer mi liquidfP^^^. Recently, Mach shocks 
from jets have been studied in Ref.'SSElEO! 

It is well knowrl^^ that a point-like perturbation (a small body, a hadron or 
parton etc.) moving with supersonic speed in the spatially homogeneous ideal fluid 
produces the so-called Mach region of the perturbed matter. In the fluid rest frame 
(FRF) the Mach region has a conical shape with an opening angle with respect to 
the direction of particle propagation given by the expressior|^ 6m = sm~^ (^) , 
where Cg denotes the sound velocity of the unperturbed (upstream) fluid and v 
is the particle velocity with respect to the fluid. In the FRF, trajectories of fluid 
elements (perpendicular to the surface of the Mach cone) are inclined at the angle 
A9 = 7r/2 — with respect to v . Strictly speaking, the above formula is applicable 
only for weak, sound-like perturbations. It is certainly not valid for space-time 
regions close to a leading particle. Nevertheless, we shall use this simple expression 
for a qualitative analysis of flow effectJ^. Following Refs.^"*^ '^^ one can estimate 
the angle of preferential emission of secondaries associated with a fast jet in the 
QGP. Assuming the particle velocity to be u = 1 and the sound velocity to be 
Cs = l/-\/3 leads to A9 ~ 0.96. This agrees well with positions of maxima of the 
away-side two-particle distributions observed in central Au-|-Au collisions at RHIC 

''Here and below quantities in the FRF are marked by tilde. 
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energies. 

4. Deformation of Mach shocks due to radial flow 

Assuming tliat tlie away-side jet propagates with velocity v parallel to the matter 
flow velocity u and u does not change with space and time, one sees that after 
performing the Lorentz boost to the FRF, a weak Mach shock has a conical shape 
with the axis along v . In this reference frame, the shock front angle 9m is again 
given hy 6m — sin^^ (if ) • Transformation from the FRF to the center of mass 
frame (CMF) shows that the Mach region remains conical, but the Mach angle 
becomes smaller in the CMF, tan6'M = J- tan 0m , where 7„ = (1 — is the 

Lorentz factor corresponding to the flow velocity u . Using the above Eqns. leads 
to the expression for the Mach angle in the CMF 

^?M=tan-i(^c.^/l^ ' 

where 

v=^^, (2) 

and the upper (lower) sign corresponds to the jet's motion in (or opposite to) the 
direction of collective flow. For ultrarelativistic jets {v 1) one can take u ~ 1 
which leads to a simpler expression 



Csls 
lu 




(3) 



with 7^ = (1 — Cj)"-^/^ . According to Eq. ([3]), in the ultrarelativistic hmit 9m does 
not depend on the direction of flow with respect to the jet. The Mach cone becomes 
more narrow as compared to jet propagation in static matter. This narrowing effect 
has a purely relativistic origin. Indeed, the difference between 9m from Eq. ([3]) and 
the Mach angle in absence of flow (lim 9m = sin~^ Cg) is of second order in the 

collective velocity u . 

The case of a jet propagating at nonzero angle with respect to the flow velocity 
is more complicated. Mach shocks become nonconical for non-collincar flows. For 
simplicity, we study only the case when the jet and flow velocities arc orthogonal to 
each other, v ^- u. Let axes OX and OY be directed along u and v , respectively. 
We flrst make the transition to the FRF by performing a Lorentz boost along the 
OX axis which leads to a jet velocity v. 

Assume a jet propagating along the path OA = vt during the time interval t in 
the FRF. At the same time, the wave front from a point-like perturbation (created 
at point O) reaches a spherical surface with radius OB ~ OC = Cgt. Two tangent 
lines AB and AC show the boundaries of the Mach regior^l with the symmetry 

''Such region exists only \lv>Ca. This condition is fulfilled li v > Cs or u > Cs hold. 
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Fig. 7. Angles of Mach region created by a jet moving transversely (solid and dashed curves) and 
coUinearly (dashed— dotted line) to the fluid velocity u in the CMF. All curves correspond to the 
case = 1/3 . The arrow marks the value u = Cs [from Satarov et al.pH 



axis OA . This region has a conical shape with opening angles 9 determined by the 
expressions sin 9 = ^ = & ~ . 

Performing inverse transformation from FRF to CMF, it is easy to show that the 
Mach region is modified in two ways. First, it is no longer symmetrical with respect 
to the jet trajectory in the CMF. The boundaries of the Mach wave have different 
angles, 9^ 9^, with respect to v in this reference frame. One can interpret this 
effect as a consequence of transverse flow which acts like a "wind" deforming the 
Mach cone along the direction OX. On the other hand, the angles of the Mach front 
with respect to the beam (OZ) axis are not changed under the transformation to 
the CMF. We conclude that, due to effects of transverse flow, the Mach region in 
the CMF should have a shape of a deformed cone with an elliptic base. Figure [7] 
shows the numerical values of the Mach angles for an ultrarelativistic jet moving 
through the QGP transversely or coUinearly to its flow velocity. 

We point out a much stronger sensitivity of the Mach angles 9± to the transverse 
flow velocity as compared with the collinear flow. 

To discuss possible observable effects, in Fig. [8] we schematically show events 
with different di-jet axes AiBi (i = 1, 2, 3) with respect to the center of a firebalE 
In the 2 — 2' event, the away-side jet '2' propagates along the diameter A2B2 , i.e. 
coUinearly with respect to the collective flow. In the two other cases, the di-jet 
axes are oriented along the chords, AiBi and A^B^ , respectively. In such events, 
the fluid velocity has both transverse and collinear components with respect to 



For simplicity we consider the case when both trigger (i') and away— side (i) jets have zero 
pseudorapidities in the CMF. 



1, 2008 16:46 WSPC/INSTRUCTION FILE paperl2 



Hydrodynamic Flow and Mach Cones 11 




Fig. 8. Schematic picture of Mach shocks from jets 1, 2, 3 propagating through the fireball matter 
(shaded circle) created in a central heavy— ion collision. Dotted arrows represent local velocities 
of the fireball expansion. Thick downward arrows show associated trigger jets. The Mach shock 
boundaries are shown by solid lines. Short-da shed lines give the positions of the shock fronts in 
the case of a static fireball [from Satarov et al.jSll, 



the jet axis. We also show how the Mach fronts will be deformed in an expanding 
matter. It is easy to see that the radial expansion of the fireball should cause a 
broadening of the sideward peaks in the A(/)-distributions of associated hadrons. 
Due to the radial expansion, the peaks will acquire an additional width of the order 
oi < 9+ — 9- > . Here 9± are local values of the Mach angles in individual events. 
The angular brackets mean averaging over the jet trajectory in a given event as 
well as over all events with different positions of di-jet axes. Assuming that particle 
emission is perpendicular to the surface of Mach cone and taking < u >~ 0.4, Cg — 
1 1 \/3 , we estimate the angular spread of emitted hadrons in the range 30° — 50° . 
This is comparable with the half distance between the away-si de peaks of the A0 
distribution observed by the STAR and PHENIX coUaborationJ^HEIIMlM], Qn the 
basis of this analysis we conclude that in individual events the sideward maxima 
should be asymmetric and more narrow than in an ensemble of different events. 
Due to a stronger absorption of particles emitted from the inner part of the shock 
(events 1, 3 in Fig. [8|), the two peaks may have different amplitudes. We think that 
these effects can be observed by measuring three-particle correlations. 

There is one more reason for broadening of the A^-distributions which one 
should keep in mind when comparing with experimental data: due to the momentum 
spread of the initial parton distributions, Ap» < 1 GeV, the di-jet system has a 
nonzero total momentum with respect to the global CMF. As a consequence, the 
angle 6** between the trigger- and the away-side jet is generally speaking not equal 
to TT , as was assumed above. Taking typical momenta of initial partons as po 7 with 
Po > 4 — 6 GeV"^*-* '^'^ '-'■^ ^"^j we estimate the angular spread as [tt — 0* | lS.p,,/pa < 
0.1 . Therefore, the considered broadening should be much less than the typical 
shift of the Mach angles due to the collective flow. 



1, 2008 16:46 WSPC/INSTRUCTION FILE paperl2 



12 Barbara Betz et al. 



0.3 r 




Fig. 9. STAR data on near-side and away— side jet correlation compared to the HSD model 
for p+p and central Au+Au collisions at midrapidity for Pxi^Trig) = 4...6GeV/c and 
PT = 2GeY/c...pT{NTrig) [from Gassing et ai.fSSESl. 



5. Angular Correlations of Jets Can jets fake the large t)2-values 
observed? 

Figure [HI shows the angular correlation of high-pT particles (pT(NTrig) = 
4 ... 6 GeV/c, PT = 2 GeV . . .pT(NTrig), \y\ < 0.7) for the 5% most central Au+Au 
collisions at ^/s = 200 GeV (solid line) as well a,s p + p reactions (dashed line) from 
the HSD-modeiS2i in comparison to the data from STAR for p+p colhsionJ^. Gat- 
ing on high-pT hadrons (in the vacuum) yields 'near-side' correlations in Au-|-Au 
collisions close to the 'near-side' correlations observed for jet fragmen tation in the 
vacuum (p+p) . This is in agreement with the experimental ob servatior l^^ l ^'^ i How- 
ever, for the away-side jet correlations, the authors of Ref.^^ get only a ^--^50% re- 
duction, similar to HIJING, which has only parton quenching and neglects hadron 
rescattering. Clearly, the observecP complete disappearance of the away-side jet 
(see Fig. [T0|) cannot be explained in the HSD-(pre-)hadronic cascade even with a 
small formation time of 0.8fm/c. Hence, the correlation data provide another clear 
proof for the existence of the bulk plasma. 

Although (pre-)hadronic final-state interactions yield a sizable (< 50%) contri- 
bution to the high-pT suppression effects observed in Au+Au collisions at RHIC, 
~ 50% of the jet suppression originates from interactions in the plasma phase. The 
elliptic flow, V2, for high-transverse momentum particles is underestimated by at 
least a factor of 3 in the HSD transport calculations*^*^. The experimentally ob- 
served proton excess over pions at transverse momenta px > 2.5 GeV/c cannot be 
explained within the CGG approach""-"; in fact, the proton yield at high-p^ > 5 
GeV/c is a factor 5-10 too small. We point out that this also holds for partonic jet- 
quenching models. Further experimental data on the suppression of high- momentum 
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Fig. 10. High pt correlations: in-plane vs. out-of-plane correlations of the probe (jet-|-secondary 
jet fragments) with the bulk {v2 of the plasma at px > 2 GeV/c) prove the existence of the initial 
plasma state (STAR-coUaboration, preliminary). 



hadrons from d-t-Au and Au-f Au collisions, down to ^/s — 20 GeV, are desperately 
needed to separate initial-state Cronin effects from final-state attenuation and to 
disentangle the role of partons in the colored parton plasma from those of interact- 
ing pre-hadrons in the hot and dense fireball. 

Can the attenuation of jets of pr > 5GeV/c actually fake the observed V2- 
values at Pt ~ 2 GeV/c? This question comes about since due to fragmentation and 
rescattering a lot of momentum-degraded hadrons will propagate in the hemisphere 
defined by the jets. However, their momentum dispersion perpendicular to the jet 
direction is so large that it could indeed fake a collective flow that is interpreted as 
coming from the early high-pressure plasma phase. 

On first sight, Fig. [TOl shows that this could indeed be the case: the in-plane V2 
correlations are aligned with the jet axis, the away-side bump, usually attributed to 
collective V2 flow (dashed line) , could well be rather due to the stopped, fragmented 
and rescattered away-side jet! However, this argument is falsified by the out-of- 
plane correlations (circles in Fig. [T0|) . The near-side jet is clearly visible in the 
valley of the collective flow V2 distribution. Note that V2 peaks at ip — tt/2 relative 
to the jet axis! The away-side jet, on the other hand, has completely vanished in 
the out-of-plane distribution (cf. Fig. [TT|). 

Where are all the jet fragments gone and why is there no trace left? Even if 
the away-side jet fragments completely and the fragments get stuck in the plasma, 
leftovers should be detected at momenta below 2GeV/c. Hadronic models as well 
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Fig. 11. Illustration of jets traveling through the late hadronic stage of the reaction. Only lets 
from the region close to the initial surface can propagate and fragment in the vacuurri ^ 1 1 62 1 65 ] 
The other jets will interact with the bulk, resulting in wakes with bow waves travelling transversely 
to the jet axis. 

as parton cascades will have a hard time to get a quantitative agreement with these 
exciting data. 

We propose future correlation measurements which can yield spectroscopic in- 
formation on the plasma: 

(1) If the plasma is a color-electric plasmgl^IMl^ experiments will - in spite of strong 
plasma damping - be able to search for wake-riding potential effects. The wake 
of the leading jet particle can trap comoving companions that move through 
the plasma in the wake pocket with the same speed (pt /m) as the leading par- 
ticle. This can be particular stable for charmed jets due to the deadcone effect 
as proposed by Kharzeev et aP3, which will guarantee little energy loss, i.e. 
constant velocity of the leading D-meson. The leading D-meson will practically 
have very little momentum degradation in the plasma and therefore the wake 
potential following th e D will be able to capture the equal speed companion, 



(2) One may measure the sound velocity of the expanding plasma by the emission 
pattern of the plasma particles travelling sideways with respect to the jet axis: 
The dispersive wave generated by the wake of the jet in the plasma yields 
preferential emission to an angle (relative to the jet axis) which is given by the 
ratio of the leading jet particles' velocity, devided by the sound velocity in the 
hot dense plasma rest frame. The speed of sound for a non-interacting gas of 
relativistic massless plasma particles is w ~ 57% c, while for a plasma 
with strong vector interactions, Cg — c. Hence, the emission angle measurement 
can yield information of the interactions in the plasma. This point will be 
discussed in the following. 
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Fig. 12. Contour plot of the laboratory energy density at t = 12.8 fm/c for an ideal gas EoS (left) 
and for a liadron gas with first-order phase transition to QGP (right). 



6. A hydro dynamical study of jet energy loss 

As men tioned in Sec. [51 the STAR and PHENIX collaborations published the 
observatiorl^SISlISn! that the away-side jet in Au+Au collisions for high-p-r particles 
(4 < pT(trigger)< 6 GeV/c, pr(assoc)> 2 GeV/c) with pseudo-rapidity \y\ < 0.7 
is suppressed as compared to the away-side jet in p-fp coUisions (see Fig. [TU)). 

This is com monly interpreted as parton energy loss, the so-called jet 
quenchinJ^nZni, One part of the back-to-back jet created in the collision escapes 
(near-side jet), the other one (away-side jet) deposits a large fraction of its energy 
into the dense matter. 

We use (3-1-1) dimensional ideal hydrodynamics, employ the (3-)-l)dimensional 
SHASTA (SHarp And Smooth Transport Algorithm)'^, and follow the time evolu- 
tion of a fake jet that deposits its energy and momentum completely during a very 
short time in a 2 fm'^ spatial volume of a spherically symmetric expanding system. 

The medium has an initial radius of 5 fm, an initial energy density of gq = 
1.68 GeV/fm and an initial profile velocity increasing by radius as v{r) = 0.02 r/R. 

The initial energy density of the jet is increased by Ae — 5 GeV/fm as compared 
to the medium and the jet material has an initial velocity of Vx = 0-96 c. We display 
the contour plots of the jet evolution at late state t = 12.8 fm/c for two different 
cases: Figure fT2l (l.h.s) shows the evolution for an ultrarelativistic ideal gas EoS. 
Here, the jet is initially located in the region between — 5fm < x < 3fm, |y| < 
0.5 fm, |z| < 0.5 fm. Figure [T2l (r.h.s) depicts the evolution for a hadron gas with a 
first-order phase transition to QGP and a jet that is initially located in the retion 
between —3 fm < a; < — 1 fm, 2.5 fm < y < 3.5 fm, |z| < 0.5 fm. 

The jet-induced shock front and a deflection of the jet in case of an EoS with 
phase transition to QGP is clearly visible. 
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7. Summary 

The NA49 collaboration has observed the collapse of both, vi- and V2-collective 
flow of protons, in Pb+Pb collisions at 40 A GeV, which presents evidence for a 
first-order phase transition in baryon-rich dense matter. It will be possible to study 
the nature of this transition and the properties of the expected chirally restored 
and deconfined phase both at the HiMu/low energy and at the forward fragmen- 
tation region at RHIC, with upgraded and/or second generatio n de tectors, and at 
the future GSI facility FAIR. According to lattice QCD resultA^, the first-order 
phase transition occurs for chemical potentials above 400 MeV. 

Ref.H9! 

shows that 

the elliptic flow clearly distinguishes between a first-order phase transit ion and a 
crossover. Thus, the observed collapse of flow, as predicted in Ref.™ 

is a clear 

signal for a flrst-order phase transition at the highest baryon densities. 

A critical discussion of the use of collective flow as a barometer for the EoS of 
hot dense matter at RHIC showed that hadronic rescattering models can explain 
< 30% of the observed flow, V2^ for pT > 2 GeV/c. We interpret this as evidence 
for the production of superdense matter at RHIC with initial pressure way above 
hadronic pressure, p > 1 GeV/bn^. 

The fluctuations in the flow, vi and V2, should be measured. Ideal hydrodynam- 
ics predicts that they are larger than 50 % due to initial state fluctuations. The 
QGP coefficient of viscosity may be determined experimentally from the fiuctua- 
tions observed. 

We propose upgrades and second-generation experiments at RHIC, which 
inspect the first-order phase transition in the fragmentation region, i.e., at 
Hb ~ 400 MeV {^/s = 4 - 12 A GeV or y w 4 - 5 at fuU energy), where the 
collapse of the proton flow analogous to the 40 A GeV data should be seen. 

The study of jet-wake-riding potentials and bow shocks caused by jets in the 
QGP formed at RHIC can give further clues on the EoS and transport coefficients 
of the QGP. 
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